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1. Introduction 

The rapidly advancing field of quantum information requires careful characterization 
of quantum states. While the state of an individual object cannot be determined, 
measurements on an ensemble can be used to reconstruct the state. It is not sufficient 
to perform ensemble measurements on a single observable, as relative phases are 
unobtainable. Rather, one requires a set of ensemble measurements on a sufficient 
number of non-commuting variables. It is then possible to transform the resulting set 
of probabilities to obtain a complete description of the state; such a procedure is known 
as quantum tomography. While the density matrix is the standard specification of a 
quantum state, it is often useful to represent states as quasiprobability distributions 
in phase space. The Wigner function W{q,p) [1112] has many useful properties as a 
phase space distribution; in particular, its marginals give the correct probability density 
of commuting observables and can be used for tomographic reconstruction of the full 
Wigner function. In deference to the Heisenberg uncertainty principle, the Wigner 
function is not a true joint probability density for position and momentum, and can be 
negative. The structure shows the non-classical nature of the state; interference effects 
apparent in pure states are smeared out in a mixed state. It is widely used to describe 
quantum states in quantum optics [3 El and has been determined tomographically in 
a number of systems, for example helium atoms passing through double slits [Hj, the 
electromagnetic field in a microwave cavity ^ and a Josephson junction coupled to a 
microwave cavity 

Expectation values take a form reminiscent of classical statistical mechanics: they 
are averages with respect to a distribution W{q,p) in phase space F: 

if) = TTifp) = j^dqdp fiq,p)W{q,p). (1) 

Here / is a function that depends only on the operator / of interest, and the Wigner 
function W is a normalized distribution that depends only on the density matrix p 
jH]. Although quasiprobability distributions W{q,p) satisfying equation (P) indeed exist 
|21in|, the only pure states for which the Wigner function is positive definite are Gaussian. 
A Gaussian smoothing over a phase space area of results in the positive-definite 
Husimi (or Q) function. However, the integral of the Wigner function over a bounded 
region in phase space may still be negative even if its area is larger than the uncertainty 
bound |?t (ini. Hence the approach to the classical limit, as a delta function on the 
classical torus, is highly singular [TT] . 

Path integration ^2] allows evaluation of expectation values in a formally similar 
way to equation ((11): 

(/) = Tr(//)) = [vqDp f[q,p]e~'^'^''^K (2) 

The integration is now over paths (g(r),p(r)) in phase space. The distribution of paths 
Q-s[q,p] jg j^Q^ general complex, as the action S[q,p] carries a geometric phase factor. 
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The purpose of this work is to relate these two formulations of a distribution, and 
show how the phases in the path integral contribute to the non-positivity of the Wigner 
function. To do this we could associate a subset of the paths to each point in phase space. 
This recalls the path-centroid method: integration of all non-zero frequency modes out 
of a configuration-space path integral maps the system to an effective classical system, 
usually determined variationally [121 CHI El. Paths x(r),0 < r < /? are classified 
according to their path centroid x = ■|Jg^a;(r)dr. The path integral with action 
S[x] reduces to an ordinary integral over c- number variables with classical effective 
Hamiltonian Hef[{x). Excursions of the path from x provide quantum corrections to 
the potential in the effective Hamiltonian. The function W{x) = exp{—f3Hcs{x)), 
regarded as a Boltzmann distribution of the variables x, forms the basis of a classical 
statistical mechanics in phase space ^H] • This distribution is an intuitive interpretation 
of the path centroid distribution if the observables represented by x are compatible, as 
is indeed the case for configuration-space path integrals, and has successfully dealt with 
quantum corrections to the dynamics of solids. However, this is not the case when the 
variables are incompatible: the action in the path integral is complex, and the resulting 
phase space distribution is not positive definite. 

The present author has previously shown a correspondence between the path- 
centroid distribution and the Wigner function W{S) for components of a free spin s 
[inillZI. This is a singular distribution, with derivatives of delta functions supported on 
spheres of quantized radius. The distribution of the path centroids of spherical polygons 
of L vertices on the Bloch sphere converges as L — oo to the Wigner function, where the 
weight of a polygon of area Q includes a Berry phase sQ. (We consider the Bloch sphere 
as embedded in M^.) For a given value of the path centroid inside the Bloch sphere, 
the dominant path is a small circle enclosing non-zero area and therefore contributing a 
phase to the Wigner function. More generally, the Wigner function of a set of operators 
can be computed as a histogram of their path averages: 

Wix) = {5Nix - x)) = lim {6Nix - P))l. (3) 

L— >oo 

Here x = {xi, . . . ,xn) are operators commuting with the Hamiltonian, a; G are 
c-number variables and the P functions of the operators are defined as distributions 
over the coherent state manifold F (e. g., the Bloch sphere S"^): 

Xi = l^dfi{^)Pi{^)\-f){-f\,t = l...N. (4) 

Pi= f Pi{l{T))dT (5) 

Jo 

is the path centroid, the time average of the P function. The first set of angle brackets 
in equation Q represents a thermal expectation value, and the second an average over 
L-vertex polygonal paths in the coherent-state path integral; is the symmetrically- 
ordered delta function in (see equation (fTHjl below). 

The above result is only applicable to compact phase space, where the excursions 
of the path are bounded for a free particle. The sign oscillations of the Wigner function 
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are a consequence of the geometrical phase arising from the curvature of the phase space 
(or coherent-state manifold). In the present work we obtain an equivalent result for fiat 
phase space, with a Hamiltonian confining the dynamics to a region of phase space. The 
Wigner function is no longer the distribution of the path centroid, but is linked to a 
distribution of the midpoint of the ends of the path. The paths are confined to regions 
where the Hamiltonian is small, and the sign oscillations of the Wigner function are 
related to the symplectic area enclosed by the path. 

In the next section we review the relevant properties of coherent states and phase 
space distributions. Section 01 contains the main result, a path integral expression for the 
Wigner function (jHSHSHI)- We apply this to obtain the Wigner function for a number 
state in Section 0J By saddle-point evaluation of the path integral we recover the 
usual WKB approximation to the Wigner function. Section El contains some concluding 
remarks. 



2. Phase space distribution 

2.1. Coherent states 

We consider a two-dimensional phase space F = with a single coordinate q and 
conjugate momentum p. The Hilbert space Ti is spanned by the number states 
{|n), n = 0, 1, . . .} of a harmonic oscillator of (arbitrary) mass m and natural frequency 
uj. The raising and lowering operators 



7 I 



a) = \ — —q p and a = \ — —q H p (6) 

satisfy the canonical commutation relations 

[a, a)] = 1. (7) 

Coherent states \a) G 7Y, with complex label a G C provide a continuous basis 
for Hilbert space: 

|«) = e°^'-"*"|0) = e-l"l'/2 J2 -^In). (8) 



n=0 ^ 

(In this paper, kets with Roman or numerical indices will always refer to harmonic 
oscillator basis states, and kets with Greek indices will always refer to coherent states.) 
Coherent states are normalized eigenkets of the lowering operator, d\a) = a\a), but are 
not orthogonal: 

(/5|7) = e-5l/^l'-|H'+/3*7. (9) 
They form an overcomplete basis for Hilbert space: 

- / d^a\a){a\ = 1, (10) 



Phase-space path integral for the Wigner function 5 
where f d^a = f°° dRea f°° dima. We note the one-to-one correspondence between 

J J — OO J —CO ^ 

the complex plane and phase space, 



muj 1 

the ket |a) is a minimum-uncertainty state centred at the point {q,p). 

2.2. The P function 

The quantum state of a system is specified by its density matrix p, a positive Hermitian 
operator with unit trace. This can be written as a weighted average of coherent-state 
projections, 

p = j d^aP{a,a*)\a){a\. (12) 

The (Glauber-Sudarshan) P function, as implicitly defined above, (see e. g. references 
jH^JE]) is the expectation value of the normal-ordered 6 operator: 

P(a, a*) = Tr [p6{a* - a^)6{a - a)] , (13) 

where 

/i2 
^gic{a--at)gic*(a-a)^ (14) 

Clearly ()12j) the P function of a coherent state is a delta function. For number 
states it is a more singular distribution. 

2. 3. The Wigner function 

Any classical state is fully specified by a non- negative distribution Wc\{q,p) over phase 
space such that the expectation value of a function f{q,p) is 



(/(g,p))ci= / dqdp f{q,p)Wciiq,p). (15) 



r 



The expectation value of an operator A can nevertheless still be written in the above 
form (fT3j) . 

Tr{fp) = J^dqdpf{q,p)W{q,p), (16) 

where the c-number symbols f{q,p) and W{q,p) are linear functions of the operators / 
and p respectively. There is a one-dimensional family of such functions, which includes 
the P, Wigner and Q functions, corresponding to distributions for normal, symmetrical 
and anti-normal ordering of the operators fH^ IT?H 1^ . We shall compute the Wigner 
function W{q,p) 012], which we shall subsequently write as a function in the complex 
plane. W{a,a*) is a linear function of p, 

W{a,a*) = TT[p62{a-a)], (17) 
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where 62 is the symmetrically ordered 2-dimensional delta function: 



7r2 



icat — ic*a 



e-|c|V2g-icatg-ic*fi 



To evaluate the matrix elements of this operator, we use the overlap and the Baker- 
Campbell-Hausdorff identity e"^"*"^ = e~''^'^l/^e^e^ (for the case where A and B commute 
with their commutator): 

(19) 

(20) 

(21) 
(22) 



giving 



[3 



52(0! — a) 



1 



d^C ; 



ica*+ic*Q!— |cp/2 , 



2ti 

d^C ; 



, _icatg-ic*a 



|7) 



ica*+ic*a— |cP/2 — ic/3*— ic*7 



27r 



-2(Q-7)(a*-/3*) 



The delta operator is also a "displaced parity" operator [21/ , which inverts a coherent 
state about the point a: 



|2a-7). 



(23) 



The pre-factor here is a geometrical phase factor, with phase four times the area of the 
triangle Oa7 (or the area of the triangle with mid-points 0,a,7 |22ll2n])- 

The Wigner function is a Gaussian convolution of the P function of the density 
matrix, 



W{a,a*) = ^ J d2/5e-2|"-^l'p(/5,/?*). 



(24) 



The Wigner function exists and is bounded for any density matrix, even if the P 
function is unbounded. Its marginals, obtained by integrating over one variable, are 
the probability densities of the conjugate variable; the Wigner function can in turn be 
reconstructed from the ensemble of measured marginals. 



3. Path integral calculation of Wigner function 

Path integration ^21 12^ is an important technique both analytically and, increasingly, 
numerically for the evaluation of a thermal expectation value as an integral over paths. 
If these paths are sequences of points in a state space they will acquire a geometrical 
phase (even in an imaginary-time path integral). Thus the evaluation of the integral over 
paths will rely on the interference between complex weights; in the numerical context 
the difficulty in evaluating such a integral by random sampling constitutes the Monte 
Carlo sign problem ^H]- Formally we can write the path integral in a similar form to 
equation fTHjl . where the integration is over paths 7(r) in the path space V : [0, 1) T: 

Tr(/p) = / P7/WPW, (25) 
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where /[7] and ^[7] are functionals defined in terms of / and p. 

To relate tliis to our discussion of Wigner functions as distributions over 
phase space, we recall the path-centroid method for obtaining effective potentials in 
configuration-space path integrals jT21II3llIlllIS]- Here each path g(r) is associated with 
the path centroid q = Jq q{T)dT and the path integral is performed over all paths g(r) 
subject to this constraint. The resulting distribution of g, interpreted as a Boltzmann 
distribution exp(— /3KfT(g)), allows definition of an effective potential V^q. 

More generally, we can classify paths by points in phase space, C : P — F, where 
C is not necessarily the path centroid. The appropriate classification, we will see later, 
is the centre of the chord linking the end-points of the path, CYi\ = (7(0) + 7(l))/2. 
Let Va = {1 & V : CYl] = a] be the space of paths associated with a point a G F 
and assume that the path integral in Va can be performed in some approximation, 
leaving a two-dimensional integral over phase space and mapping functionals of paths 
onto functions in phase space: 

j^V^fb]= j^d'afia) (26) 

where 

/»= / P7/[7]= / V^f[^Ma-C[^]). (27) 

The contribution of the subspaces in this integral can then be interpreted as a 
quasidistribution over phase space; any non-positivity arises from interference between 
the geometrical phases of the paths, the sign relating to the distribution of areas of paths 
in V^. Of course, the integrand y4[7]p[7] in the path integral will not in general 
separate after path classification. 

The relation of the phase of the Wigner function to areas in phase space is not 
new, as Schleich has extensively reviewed jS]. Berry ^l] has established a relationship 
between the Wigner function and areas in phase space in a semiclassical context: a line 
segment is drawn that intersects the classical torus in two points equidistant from (g,p), 
and the phase is then the area enclosed by the line segment and the classical torus. 
Ozorio de Almeida [22] has obtained a path integral representation for the Wigner 
function of a product of operators, and hence for the evolution operator, in terms of the 
star product of Wigner functions. Here the geometric factor is the symplectic area of a 
polygon whose mid-points are the phase space points on the path. 

Our approach is complementary to and simpler than that in references |22], and 
relates the Wigner function to a path integral of a product of P functions. There is a 
choice of ways of evaluating a coherent-state path integral, using different symbols for 
the operators concerned The present approach starts from the usual time-slicing 
of the density matrix. Suppose we wish to find the density matrix p = exp(— /3ff) 
of a Hamiltonian whose P function P//(a,a*) exists. This class includes (but is not 
restricted to) Hamiltonians polynomial in a and ^HlCSl- The density matrix can be 
time sliced as 

p= lim (i-(3H/l) . (28) 
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(There is no additional difficulty in introducing a Trotter decomposition of the 

Hamiltonian into non-commuting terms.) Replacing the Hamiltonian by its P 
representation (|T^ and re-exponentiation defines the Lth approximant to the density 
matrix as 

hiP) = n / (29) 
1=1 ^ 

= iPiiP/L))'. (30) 

It is useful to ensure that the approximant obtained numerically from finite time-slicing 
avoids unphysical pathologies [2^]. Since each approximant is (up to normalization) a 
physically realizable density matrix, the Wigner function is the limit of a sequence of 
physical Wigner functions. 

The Wigner function, following Eqs. (|T7jl and (j221), is 

r ^ 

WL{a,a*) = Tr / l[d\62{a - a)\jL) {1l\1l~i) ■ ■ ■ (72I71) (7i|e^^^'- (31) 
^ 1=1 

= I ^dSe-''^t^•"l-^^'-''«(^"^'*^/^ (32) 
1=1 

where the geometrical action 5'L[7,a] = 5']^^*'^[7] + Sf^'^[y,a\ has two terms: the 
Bargmann invariant 5'^'^*'^ depends only on the geometrical properties of the path, 

^-sr\,] ^ (^^|^^_^) . . . (72|7i)(7i|7l), (33) 



7/ - 7/_i| + ilm7«7;_i , (34) 



and Sf^'^, connecting the end-points of the path to a, is obtained from ()22|) as 

5r'i[7,a] = 2(a-7L)(«*-7i*)- (35) 

(Equation ()34|) uses "periodic boundary conditions" 70 = 7l for notational convenience.) 
The total geometrical action can then be manipulated into the form 

2 



Re^L[7,a] = ^ ^ |7« - 7i-i|^ + 2 



2 

1=2 



7i + 7l 
a 



(36) 



Im S'l[7,q;] = Im^7z7;*_i + 2 Im(a - 7L)(a* - 7^). (37) 
1=1 

Figure^lshows the geometrical interpretation of the above action. The path is classified 
by the mean of its end-points, 

7 = C[7] = ^4^. (38) 

The real part of the geometrical action is equal to half the sum of the squares of the 
line segments shown in bold in figure ^ that is, the L — 1 "internal" segments of the 
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Figure 1. Geometrical properties of path for evaluation of Wigner function at a (for 
L — 4). The real part of the geometrical action depends on the sum of squares of the 
bold line segments; 7 is the mid-point of the line 7171, joining the ends of the path, 
and is the mid-point of the line a/3. The imaginary part is proportional to the shaded 
area, where a lies on the central line of the rectangle DEFG. 

path (excluding 7/, to 71) and the hne from a to 27 — a, the first and second terms in 
dnSl) respectively. The second term suggests a Gaussian smoothing of the distribution of 
chord midpoints: the centre of the path end-points is close to a, and the delta function 
in is effectively replaced by a Gaussian. The dynamical part of the action favours 
regions of phase space where the P function of the Hamiltonian is small. The imaginary 
part is twice the area shown shaded in the figure: the area enclosed by the path plus 
the area of the rectangle with the line segment 71 to 7^ forming one side and a on the 
line of symmetry parallel to this side. Since the time-reversed path changes the sign 
of the area but leaves the lengths invariant the Wigner function is real. Furthermore, 
for a outside the energy surface, the dominant path has all points clustered beneath a 
and encloses zero area. The Wigner function is therefore oscillatory inside the energy 
surface and monotonically decaying outside. The real part of the action also represents 
a ferromagnetic XY model on a one-dimensional chain with L spins, a local anisotropic 
potential Ph/L allowing for amplitude fluctuations, an antiferromagnetic bond between 
the ends 71, 7l of the chain, and an external field (due to a) acting on the end spins. 

4. Example: number state 

As a concrete example, we evaluate the Wigner function of a number state. To do 
this, we define a family of density matrices Pl{N) interpolating between a Poisson state 
Pi{N) of mean occupation and a number state poo{N). Here L is a positive integer 
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and is a non- negative real number. We first define pi{N) to be a phase-randomized 
coherent state of magnitude |a| = a/ZV fI7\ : 



p,(N) = j'^^ v^e^^) (v^e^^ (39) 

From equation (jHJ this is easily shown to be a Poissonian mixture of number states of 
mean occupation N: 

Jo 27r ^^^^ vV^ 

°° AT" 

= ^"''E^NH- (41) 

n=0 

The P function of the distribution is localized on a circle, 

Pi(a,a*;Ar) = 25(|a|2-A^), (42) 
and the Wigner function, obtained from the integral ()24|) . is 



Wi{a, a*;N) = -e-2|°l'-2%(4yiV|«|), (43) 

71 

where Ii, is the modified Bessel function. The P function and Wigner function are 
non-negative for this state. 

The family of density matrices is then defined by radially squeezing the above state, 
reducing the number fluctuations: 

^ n (/'^ ^) I7l)(7l|7l-i) ■ ■ ■ (72|7i)(7i| (45) 



Zl 

with normalization factor 

Z^(iV) = Tr(pi(iV))^ = e-^^5^(^) . (46) 

The coefficient of \n){n\ dominates in the sum if?2<A^<'ri + l. Thus, for any 
non-integer A^, Pl{N) converges onto the number state \n){n\ in the limit L oo, 
where n is the largest integer less than A^. We choose = n + ^ to obtain the fastest 
convergence onto the state |^^)(^^|. This corresponds to the classical energy surface at 
the energy eigenvalue, 

P\^ = nu;(n + l). (47) 



2m 2 V 2^ 

The P function of the number state is singular, while the Wigner function can be 
expressed in terms of Laguerre polynomials [U HH] : 

P^ia,a-;N) ='—g^SM), (48) 
»'o.(a,o*;iV) = i(-l)"e-2l"PL„(4|a|2). (49) 

71 
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Thus the P symbol becomes increasingly singular with L, while the Wigner function 
remains bounded but develops sign oscillations. 

It is useful to ensure that the approximants (jHSI) to the number state Wigner 
function refer to physically realizable states |2HI. This is indeed the case here, as the 
density matrices (jl^ constitute a sequence of thermal density matrices for the (non- 
polynomial) Hamiltonian 

H{N) = N + \nn\-n\nN (50) 

POO 

= N + ln e-^+^'^^^^/^^rfa; (51) 



at temperatures kT = 1/L,L = 1, . . . , oo. The path integrals to be obtained are 
therefore exact expressions for the thermal Wigner function of a physically admissible 
Hamiltonian. 

While the Hamiltonian ()5H) does not have an obvious physical interpretation, for 
large N we can use Stirling's approximation to expand about the minimum to give a 
Hamiltonian, quadratic in the number operator, that might for example describe the 
charging energy of a Cooper-pair box 

HiN) ^ - In 2nN + ^ (52) 

We now calculate the path integral as defined in section El For convenience we will 
use polar coordinates 

a = se'*. (53) 

The path vertices are restricted to the circle 

7i = re^^' (54) 

where r = The symmetry allows us to take = in the following (or to define 

the 6i relative to 0). The Wigner functions ()32|1 simplify to 

with the action defined as in equations (|H1HH7|) : 

L 

SlIi, a] = Lr^ + 2s^ -r^Yl e'^''"^"''^ + 2r^e"^'^-''^ - 2rs{e-''' + e''^). (56) 

1=1 

The real part of the action is that of a ferromagnetic XY model on an open chain of 
L spins, with an antiferromagnetic bond between the end spins and an external field 
acting on the end spins. The imaginary part is related to the area enclosed. 

The saddle-point approximation to the integral ()55|) is given in the appendix. The 
result is a Wigner function oscillatory inside the energy surface |a| = \/n + 1/2 and 
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Figure 2. The Wigner function for tlie number states |1) and |10). Tlie daslied line 
represents the exact Wigner function The bold line represents the saddle-point 

approximation derived here jSTJ; the overall amplitude has been adjusted to agree with 
Berry's result |IT|. The non- negative dotted line represents the Wigner function for 
the Poisson state pi{n+ i) 



decaying outside: 



(2n+l)cos-i|| l"l ^ j-2|a|^n+§-|a|2-7r/4 



W^{a,a*;n+-) oc { 



1/4 



\a\ < \/ n + ^ 



(57) 



exp 



(2ra+l) cosh" 



1 kl 



= -2|a|^|a|2-„-i 



\a\ > \ n + 



The functional form is identical to the Wigner function evaluated from the WKB 
approximation to the wave function of \n) [TT] . 



W 



WKB/ 



COS ( {2n + 1) cos 



V^2 



2\a\\/n + 



\a\ 



- 7r/4 



(58) 

(7r3/2)i/2(l«l'(^ + l/2- |a|2)'^" 
although the saddle-point approximation does not yield the correct numerical coefficient. 

Figure El shows the Wigner function and its saddle-point approximation; the only 
adjustable parameter is the overall amplitude, which is taken from equation (j3H)) . The 
fit is very good apart from singularities at the origin and on the energy surface. The 
fit to the Wigner function for \n) requires fixing the parameter N = n + ^. An exact 
evaluation of the path integral for WL{ci,a*; N) should converge to the same Wigner 
function for n < N < n + 1. The exact Wigner function Woo{c(, a*; A^), must have 
discontinuities at every integer value of A^. The saddle-point approximation smears 
out these discontinuities, so that the Wigner function continuously interpolates between 
that for the number states |n — 1) and \n). 
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5. Conclusions 



The path integral formalism developed here has two pm'poses. Firstly, it provides a 
possible method for the calculation of the Wigner function for a wide range of density 
matrices in phase space. This includes, but is not restricted to, thermal distributions 
resulting from a polynomial Hamiltonian; a simple example might be a Cooper-pair 
box. The example of a number state shows that good results can be obtained in a 
saddle-point approximation. 

Secondly, it provides a further geometric interpretation of the sign oscillations of 
the Wigner function in terms of the area enclosed by the dominant paths associated 
with a point {q,p) in phase space. To find the quasiprobability density that a particle 
is at {q,p), we seek open paths in phase space associated with this point. The action 
has three terms, confining the path to low-energy regions of phase space (more strictly, 
where the P function of the Hamiltonian is small), minimizing the length of the path 
and placing the mid-point of the ends of the path close to the point For points 

outside the energy surface the dominant path will enclose no area, and the Wigner 
function will be positive; for points inside the energy surface the path encloses a finite 
area, and the sign of Wigner function will depend on the area enclosed. 



Appendix: Saddle point evaluation of the path integral 

Here we outline the evaluation of the path integral for the Wigner function (j37|) in the 
saddle-point approximation. The action 5'l[7,q;] is an analytic function of {9i} 
(but not of {7/}), 

SLb, a] = (a) + IY1 S^l\c^)UOi - ^^Dl^-^ " + • • • , (A-l) 



2 

Im 



where the action is stationary, dSi/dOi = 0, at the saddle point 61 = 6f\ The derivatives 
are 



i [_^2 ^^Ke,-9,) ^ + 2rse-''^] , (A.2) 

^^m-i-eo _ ei{e.-e.+i)) , 1 < / < L, (A.3) 

i \r^ (e'(^^--^^) + e^(^^-^^)) - 2rse'^-l . (A.4) 

Setting the derivative ()A.3|) to zero shows that the angles are equally spaced around an 
arc of angle 26, 

9,'°' = ( ^-^^^) 0, (A.5) 



L-l 

where substitution in equation ()A.2j) or ()A.4jl gives an implicit equation for 0: 
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The saddle point has a non-zero imaginary part for finite L, and for s > r the saddle- 
point value of 6 is purely imaginary. However, as L ^ oo the saddle point becomes 

s = r cos 9, (A. 7) 

so that the chord joining the end-points of the arc passes through the point a. 

The action at the saddle point follows from substitution of Eqs. flA.SHA.G]) into 
equation 

^2 



Sf\a) = Lr' (l_e-2^''/(^-i)) + ^(e- 
= ir^ (2^-sin2^) + 0(L"^) 



2i6»(L+l)/(L-l) 



) 



(A.8) 
(A.9) 



2 / 1 S 

= 2ir cos - 
V r 


— sV — 


) + 0{L- 








(A.IO) 


The second derivative at the saddle point is ^ 


liven by the L 


X L matrix (for 


L > 2) 






■■■ 





t \ 








-1 2 


_1 ... 












(^)^^2^-2i./(L-l) 


-1 


2 ■ • ■ 










(A.ll) 







■■■ 


2 


-1 








^ t 


■•■ 


-1 


2 ) 







where 



t 



^2iLe/{L~l) 



The determinant of the matrix is straightforward to evaluate: 



det4^^(a) 



t-^ [{1 + L) + 2t + {1 - L)t^] . 



The Wigner function in the saddle-point approximation becomes 

(det5f(«))-^/^ 

^(l-e^'' + 0(L-i))- 



e-4°'(") 
(27r)V2^^(r2 

gir2(sin 26l-26»+0(L"l)) 
LV2(27r)V2^^(^2)^ 



-l/2gie_ 



Substituting for 6^ from equation ()A.7|) gives 

p2isv^r^— s^— 2ir^ cos"-'- (s/r)+i7r/4 

Wl^ia, a*; r') ^ 



1/2 ■ 



(A.12) 
(A.13) 

(A. 14) 
(A.15) 

(A.16) 



r^-i2LV2(27r)V2Zi(r2) (s^r^ - s^) 
Adding the time-reversed saddle-point 6 —9 and setting = 72 + 1/2 gives our result 



()57|) for the number-state Wigner function 



cos {{2n + 1) cos-i (^7=t) " 2|a|^n + i - |«|2 - 7r/4 + 0(L"i)^ 



1/4 



(27r)V2Li/2^^(,^ + 1/2) (|a|2(n + i 

As expected, the area enclosed by the path largely determines the phase, with a small 
phase shift from the saddle-point expansion. 



(A.17) 
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For s > n + i (outside the energy surface) there is one imaginary saddle point, 



9 = i cos ^ - 

r 

and the Wigner function becomes 



(A.18) 



Wl^{a,a*;n+ 



^ exp (^(2n + l)cosh-i (^ft^) - 2|«|-\/|ttP-^ - | + 0(L-i)^ 



2(27r)V2^i/2^^(^ + 1) (|«|2(|«|2 

For large n we can use Stirling's approximation for n! ~ (^27r)-'^/2g-n+i/i2nj^n+i/2 ^.^ 
obtain 



.(A.19) 



^L(n + l/2) 



27rfn+- + 0(n-i; 



-L/2 



(A.20) 



consistent with the standard deviation y^H- | of the Poisson distribution ()41|) . and 
hence 



PFf (a,a*;n+ 



cos I (2n+l)cos-l I ^ j -2|a|^n.+l/2-|a|2--7r/4 



I a I - 



1/4 



lal < -i/ra + i 



.(A.21) 



exp 


(2n+l)cosh-i 




)-2|a|^|a|2-n-ij 






2i^/^(i+2ir 










1/4 



The functional form agrees with the WKB expression (j3Hj) of Berry up to a constant 
factor weakly dependent on n and L (but divergent in the L — oo limit). 
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